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DEVELOPMENT OF A MODEL FOR DISTRIBUTED SATURATION 
CHARACTERISTICS IN INTEGRATED DEVICES* 

By W. W. Happ 

National Aeronautics and Space Administration 
Cambridge, Massachusetts 
and 

L. B. Dickson 
Motorola Semiconductor Co. 

Scottsdale, Arizone 

SUMMARY 

The characteristics of an integrated circuit transistor, 
operating in the saturated mode, cannot be determined from 
simple lumped parameter relationships. Due to the geometry 
of these transistors, the series collector resistance must 
be calculated using conformal mapping techniques. 

To accurately determine the saturation characteristics, 
distributed models are developed for the five regions of the 
transistor. Using these models as a guide, difference equa- 
tions describing each region are formulated. Insight into 
the operation of the transistor can be gained by examining 
the internal current densities predicted by the model. 

INTRODUCTION 

The saturation voltage specification is often the most 
severe specification placed upon an integrated circuit 
transistor, and many times the ac as well as the dc-voltage 
performance is limited by this specification. To truly 
optimize the design of an integrated circuit transistor, a 
method for calculating the saturation voltage, V ce ( sat \, is 
necessary. (All symbols to be used are defined in the Appendix.) 

The saturation characteristics of the integrated circuit 
transistor in Figure 1 depend upon two-dimensional current 
flow. The following must be considered: 

(1) The effect on the emitter-base potential, V e , of the 
lateral current flow in the distributed base 
resistance . 

(2) The combined effects of lateral current flow in the 
base and collector regions on the collector base 
potential, V c . 

(3) The effect of the planar base and collector contacts 
on the series base resistance, Rgi , and series 
collector resistance, R sc . 

*The work reported here was published in part an an M.S. thesis 
of L. B. Dickson at Arizona State University. 



This investigation considers all of these effects. 

Previous related investigations by Ebers and Moll (ref. 1) , 

Gosh (ref. 2) , and Hauser (ref. 3) have only partially 
treated this problem. Dr. F. A. Lindholn of the University of 
Florida and Dr. J. Staudhammer of Arizona State University have 
significantly contributed ideas and valuable suggestions to this 
work. 


ASSUMPTIONS FOR THE DISTRIBUTED MODEL 

A one-d.imensional model for the intrinsic transistor has 
been postulated by Ebers and Moll (ref. 1) . This model gives 
the black-box performance of a three-layer transistor in terms 
of the junction potentials, terminal currents, and A(e,c) 
parameters, as defined in Eqs . (1) and (2) in Figure 2, using the 
terminology defined in the Appendix. 
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The intrinsic saturation voltage, Vce(sat)*/ ma Y be 
obtained from this model by solving Eqs. (1) and (2) for V e and 
V c . 


then 

^ce (sat) * 
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V ce (sat)* is Proportional to log aj . For small values of ag , 
such as those experienced in integrated circuit transistors, 
V ce (sat)* i s sensitive to changes in ai. 

The expression for V ce ( sa -|-)* neglects the effect of R sc . 
The standard approximation defining the total saturation 
voltage 
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is valid only for conventional transistors with one-dimensional 
flow. 
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OHMIC BASE CONTACT 



Figure 1.- Cross section of one half of an integrated circuit 
transistor showing the structure, distributed re- 
sistances, basic regions, and X axis 


-v E + + v c - 



E B C 


Figure 2.- Currents and voltages for the Ebers and Moll model 
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The lateral flow of collector and base currents in 

integrated circuit transistors prevents the separation of 

the terms in such a simple manner, and for such a transistor, 

V , V , V , . *, and a T are functions of R and the base 

G G Ce VS3."Cj -1- sc 

resistance, R b . In fact, V c and V e become functions of the 

lateral distance in Figure 1 and cannot be treated as constant 

across the transistor (refs. 2, 3). 

To examine the operation of the transistor under such 
conditions, a distributed model with one-dimensional transis- 
tors interconnected with appropriate resistors is developed. 

Gosh (ref. 2) developed a similar model for operation in 
the active mode, but did not apply it to saturated mode of 
operation, and, thus, did not include the effect of R . 

It is assumed that the transistors to be modeled in this 
manner will be gold diffused to improve their switching per- 
formance. This technique essentially eliminates the active 
parasitics created by the substrate-collector junction of an 
integrated circuit transistor (ref. 4). Therefore, the effect 
of this junction is simply to constrain the collector current. 

This constraint will be considered in the calculation of R„^ 

s c 

and accurate results may be obtained. 

DEVELOPMENT OF THE DISTRIBUTED MODEL 


The transistor shown in Figure 1 may be divided into 
five regions : 


Region I : 

Region II: 

Region III: 
Region IV: 

Region V: 


The intrinsic region consisting of the 
cross section which includes the emitter 
The area between the emitter and base 
contact 

The equipotential base contact region 
The region between the edge of the base 
and the base contact 

The region between the base and collector 
contact 


Regions II, III, and IV will be referred to as the 
extrinsic regions. Each section of the intrinsic region 
(Region I) , represents a classical one-dimensional transis- 
tor with the emitter area equal to the collector area. The 
sections in the extrinsic regions are collector base diodes. 
In Region III, the effective base resistance is zero because 
of the ohmic contact, and in Region IV the base resistance 
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is assumed to be zero for simplicity. The validity of this 
approximation is discussed later. 

Figure 3 shows the distributed model for all regions of 
the transistor. Using this figure as a guide, V ce ( sat j may 
be found by summing the terminal voltages. 


V 


ce (sat) 


= V + 
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where the above voltages are functions of the terminal currents 
and the "A" parameters of each region. Closer examination 
shows that the voltages in Eq. (5) are also strong functions 
of the lateral collector and base currents; therefore, V ce ( sa ^-) 
determined from Eq. (5) will be accurate for "top collector" 
transistors . 


The values of the voltages in Eq. (5) can be found by 
writing difference equations rather than differential equa- 
tions using the distributed model as a guide. For Region I 
(Figure 4), the current relationships are given by: 

Ifao “ (1 - V A ee (ex ? + (1 ' “l> *oc [ e *P ^ -l]<6> 
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V el - V ao + I bo *b (8) 

V cl - V co + r bo R b - R o < 9 > 


Subscripts 0 and 1 refer to the initial and first region 
of the model shown in Figure 4. Typically, the change in 
voltage from one section to another is in the order of few 
millivolts, but because the currents are exponential func- 
tions of these voltages, this change in voltage is important. 
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igure 4.- Distributed model of region I - the intrinsic region 






The equations for the nth section in Region I are: 
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(13) 


then, V eN , V cN , l faN , and I cN can be found by setting n = N 
The equations governing the performance in Region II 


are : 
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The equations governing Regions III and IV are identical 
to those for Region II with R^ = 0. 

With the equations written above, the necessary voltages 
to find v ce ( sat ) are available. 
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Substituting Eqs. (17), (18), and (19) into Eq. (5) 

N n 


V , , > = V - V 

ce(sat) eo co 


'1 

+ E 


„ Z ck R o + I cN, R c. (20 » 
k = 0 11 


The terminal base and collector currents are evaluated at 
n = and similar equations can be written for them. 

With the results from the above equations, the necessary 
currents and voltages are given in terms of V c q and V e g and 
the A(e,c) parameters. With a judicious choice of V c g and 
V e g and the use of a computer, V ce ( sa t)can be calculated for 
any desired terminal current. 

THE EFFECTIVE INVERSE CURRENT GAIN 

The effective aj of an integrated circuit transistor 
is a function of the lateral position of the current being 
injected from the collector into the base. Figure 1 shows 
the collector junction is much larger than the emitter 
junction, and current injected long distances away from the 
emitter will not contribute to the inverse emitter current 
and thus reduce ai. 

Assuming an inverse current gain, aig, for Region I 
and an inverse current gain of zero for the extrinsic regions. 
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then the results of the previous section can be used to cal- 
culate the effective aj. If Ifcj_ is defined as the total base 
current which originates in the collector and I c ^ as the 
total current injected out of the collector, the alpha ef- 
fective inverse is given by: 


a Ief f ~ 


(1 



where I . is given by: 
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where the summation extends over all regions of the transis- 
tors , and 
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THE SERIES COLLECTOR RESISTANCE 

The series collector resistance, R sc , may be approxi- 
mated to be the resistance between two conducting planes as 
shown in Figure 5b, with one contact representing the col- 
lector N+ diffusion and the other the emitter diffusion 
(Figure 5a). Assuming charge neutrality (i.e., V^p = 0) in 
the collector area, mapping techniques (Eqs. (5), (6), and 

(7)) may be used to find the capacitance for this geometry: 
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Figure 5.- Shows the similarity of the cross-section of the integrated circuit 
transistor (Part a) to the geometry, which was conformally mapped 
by Kaiser and Castro 





and Z in the length of the stripe perpendicular to the 
directions shown. K and K' are the complete elliptic integral 
of the first kind of its conjugate respectively; i.e. : 

n r U 2-] 1/2 

K. = K [l - (^ ] 

Utilizing the analogy between capacitance and resistance, 
RC = p ee 0 , Smythe (ref. 8) , where R and C are the lumped 
equivalent resistance and capacitance, respectively, and p and 
ee 0 are the distributed resistivity and permitivity of the 
material, these results can be used to calculate the series 
resistance for an integrated circuit transistor. The series 
collector resistance for a transistor with emitter length Z 
is : 


R 


SC 


p £ £ 



K (U 1 /U 2 ) 

K 1 ‘ U l /U 2> 


(28) 


The analysis has been two dimensional , but for most 
transistors of interest, Z is >> L,d and the end effects can 
be neglected. If Eq. (28) is used when Z = L, worst-case 
value results. 

LIMITATIONS OF THE MODEL 

The model is accurate for low current ranges, but has 
anticipated limitations at high current levels. Important 
factors not included in the model which affect the results 
at high current levels are: 

(1) At high current levels, alpha will decrease because 
of the electronic junction motion which increases 
the base width, (refs. 9, 10), and, thus, reduces 
the base transport factor and the emitter efficiency. 
For large currents, conductivity modulation also 
reduces the emitter efficiency, and the net result 

is that, at high current levels, alpha will decrease 
rapidly . 

(2) The resistivity and associated resistances for the 
base and collector regions are inversely related to 
the majority carrier concentrations in these regions. 
At high current densities, the number of majority 
carriers is radically increased to maintain charge 
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neutrality, and the resistivity and the values of the 
resistances are reduced. 

(3) A third limitation to this model results from the 
linearization of the series collector resistance. 

It was shown that this resistance is effected by the 
crowding of the current near the collector contact and 
emitter, and the method for accurately calculating the 
resistance is given especially for high currents. The value 
of collector resistance used in each section of the model was 
found by dividing the total R s by the number of sections. 

This linearization assumes a constant distributed collector 
resistance. This assumption introduces a second-order 
inaccuracy which becomes important at high current levels. 

APPLICATIONS OF THE MODEL 

The model is a distributed version of the non-linear 
model proposed by Narud and Meyer (ref. 9) and will accurately 
predict operation in the active mode. In the active mode, one 
can predict the emitter current density for a given total 
emitter current. As was shown by Kirk (ref. 10) , the maximum 
number of carriers passing through the collector-base deple- 
tion region must be kept less than the ionized impurity dis- 
tribution in that region for the transistor to operate. 

The number of carriers passing through a region as related 
to the current density may be determined. This model gives the 
current density in any region, and thereby, provides a means to 
exsure that the maximum allowable current density will not be 
exceeded . 

The model may be expanded to give the transient response 
of an integrated circuit transistor. It would be necessary to 
include both junction and diffusion capacitances for each sec- 
tion of the model. An analysis based on the model of the col- 
lector region would give the storage time of the transistor; 
however, the relationships are nonlinear and a computer solu- 
tion would be necessary. 

PRINCIPAL CONTRIBUTIONS OF MODEL 

The model here proposed, and the techniques developed to 
utilize it, offer significant aids to the device designer. 

(1) Appropriate formulas were devised for the calculation 
of the series collector resistance of integrated 
circuit transistors. 
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(2) The expressions for the saturation voltage and 
effective ct^. were developed from a distributed 
model. 

(3) Perhaps the most important result of this analysis 
is the insight into the operation of the transistor 
that has been gained by examining the potentials and 
currents inside the transistor as a function of 
lateral distance. This insight aids in the opti- 
mization of transistors for a given function. 

An experimental verification of the model .here developed 
shows excellent agreement with theory, and has been discussed 
elsewhere (refs. 11, 12). 
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APPENDIX 


Glossary of Terms 


A basic parameter used to model transistors - defined 
in Eqs. (1) and (2) 

The common base current gain of a transistor operating 

in the normal mode, i.e., a = I /I , V r = constant 

n c e c 

The common base current gain of a transistor with the 
collector acting as the emitter, i.e., in the inverse 
mode. 

Permitivity 

Diffusion coefficient of minority carriers in a 
semiconductor 

Base current 

Collector current 

Emitter current 

2 

Current density in amps/cm 

Boltzman's constant = 8.63 X 10 eu/°K 

Reference voltage at 25°C of value - 26 millivolts 

The complete elliptic integral of the first kind and 
its conjugate 

General section of distributed model, e.g. , Eq. (12) 

Diffusion length of holes in N-type semiconductors 

The section immediately under the base contact 

The section at the edge of the emitter 

The section at the edge of the base-collector 
junction 

The equilitrium concentration of minority carriers in 
the base of a PNP transistor 



B 


'B- 


Rr 


P 

V 


SC 


ce(sat) 


V ce(sat) 

V 

e 

V 

c 

w 


-19 

Electronic charge - 1.6 X 10 es u 

Distributed base resistance under the emitter (in the 
active region) 

Distributed base resistance in the enternal region 

The lumped resistance external to the collector-base 
junction 

The total series collector resistance 
Resistivity 

*The saturation voltage of an ideal one dimensional 
transistor with no series resistance 

The saturation voltage 

Emitter-base potential at any point 

Forward biased collector-base potential at any point. 
The base width of the transistor in the active region 
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